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Abstract. 

From small-angle neutron scattering studies of the flux line lattice (FLL) in 
CeCoIns, with magnetic field applied parallel to the crystal c-axis, we obtain the 
field- and temperature-dependence of the FLL form factor, which is a measure 
of the spatial variation of the field in the mixed state. We extend our earlier 
work [A.D. Bianchi et al. 2008 Science 319 177] to temperatures up to 1250 mK. 
Over the entire temperature range, paramagnetism in the flux line cores results 
in an increase of the form factor with field. Near Hc2 the form factor decreases 
again, and our results indicate that this fall-off extends outside the proposed 
FFLO region. Instead, we attribute the decrease to a paramagnetic suppression 
of Cooper pairing. At higher temperatures, a gradual crossover towards more 
conventional mixed state behavior is observed. 



PACS numbers: 71.27.-|-a,74.25.Qt,61.05.: 

1. Introduction 



The heavy- fermion superconductor CeCoIns continues to excite great interest, because 
it shows strong Pauli paramagnetic effects [Hll] and also the close proximity of 
superconductivity to a quantum critical point [3ll4] . It has a superconducting transition 
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temperature, ^ 2.3 K in zero field [5], with a d^2_y2 order-parameter [BJ[7]. At low 
temperature, the transition to the normal state is first-order [HH], showing that the 
superconductivity is suppressed by coupling of the field to the anti-parallel spins of 
the singlet Cooper pair (the Pauli effect) rather than the more common coupling 
to the orbital motion of Cooper pairs in the mixed state (the orbital effect) [8l|9]. 
The dominance of the Pauli effect, combined with the quasi-2D structure and super- 
clean crystal properties {£ ~ 1000^ [lOj). means that the stringent requirements for the 
stabilisation of the inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state are 
satisfied [TIHn] . Numerous studies report experimental signatures that are compatible 
with the formation of FFLO [T4HI8] , though an unambiguous microscopic observation 
remains elusive. Recently, both nuclear magnetic resonance (NMR) [19] and neutron 
diffraction studies |2Q] have provided microscopic evidence for the stabilisation of 
field-induced long range antiferromagnetic order for fields applied within the basal 
plane. This magnetically ordered phase (termed 'Q-phase' in [50]) exists only within 
the superconducting mixed state, and occupies the same high field, low temperature 
region of the superconducting (H, T) phase diagram as proposed for FFLO. It remains 
unclear whether this ordered phase replaces, or coexists with, a non-standard FFLO 
state. 

Small-angle neutron scattering (SANS) studies of the superconducting mixed 
state of CeCoIn5 have determined the flux line lattice (FLL) structure and orientation 
for the major part of the superconducting phase diagram with magnetic field applied 
parallel to the crystal c-axis (H || c) [2TH24] . In addition, SANS can be used to 
determine the FLL form factor (FF) which is a measure of the field contrast in the 
mixed state, i.e. the difference between the local magnetic induction in the fiux line 
cores and at the field- minima between them [^T^El] . Most notably, an anomalous field- 
driven enhancement of the FF was observed at 50 mK upon approaching Hc2 (T = 0) = 
4.95 T [23]. We proposed that the increasing FF arises from paramagnetic moments 
induced in the flux line cores where the antiparallel alignment of electron spins in 
Cooper pairs is suppressed, a phenomenon that is consistent with the predictions of 
recent numerical calculations [IS]. Just before entering the normal state, an abrupt 
fall of the FF was observed, which we speculated was associated with the formation of 
the FFLO state [23]. However, the numerical study [25] reproduced a similar reduction 
of the FF at high field, apparently unrelated to the onset of the FFLO state. Hence, 
the physical origin of the decreasing FF just below Hc2 remains unclear. To address 
this question, we have used SANS to investigate the field-dependence of the FLL FF 
in CeCoIns for temperatures up to 1250 mK. For all temperatures in this range, the 
FF initially rises with increasing field, reaching a maximum before it starts to fall 
prior to entering the normal state. Our results show that while we are unable to rule 
out a contribution to the field-dependence of the FF due to a FFLO-type state, the 
dominant contribution to the high field fall in the FF at low temperature is due to a 
paramagnetic suppression of the Cooper pairing. 

2. Experimental 

Our experiments were performed on the SANS-I instrument at the Paul Scherrer 
Institut. The samples were mosaics of co-aligned CeCoIns plate- like single crystals 
(with small dimension along the c-axis) grown from an excess indium flux [5]. Most 
results reported here were obtained using a sample of total volume 0.2 x 11 x 14 mm'^ 
and of total mass 250 mg. Thin samples were necessary due to the strong neutron 
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absorption of indium. The sample was mounted in a dilution refrigerator with base 
temperature ^ 50 mK which was inserted into a horizontal field cryomagnet, and 
oriented so that the c-axis was < 2° from the field direction. Neutrons of wavelengths 
A„ = 5 — 6 A were selected with a wavelength spread AA„/A„ = 10%, and were 
incident approximately parallel to the field; the diffracted neutrons were collected by 
a position-sensitive multidetector. At all fields the FLL was prepared by field-cooling 
through Tc. Diffraction measurements were performed as the sample and cryomagnet 
were rotated together to carry the FLL diffraction spots through the Bragg condition. 
Background subtraction was performed using measurements taken either above Tc or 
after zero-field cooling, with no noticeable dependence on which method was used. 

3. Results and Discussion 

The local field in the mixed state may be expressed as a sum over its spatial Fourier 
components with indices (/i, fc), and scattering vectors qi^i^ belonging to the two- 
dimensional FLL reciprocal lattice. The form factor (FF) at wavevector q^^^ is the 
magnitude of the Fourier component F(q^^). The value of the FF is obtained from 
the integrated intensity of a FLL Bragg reflection as the FLL is rotated through the 
diffraction condition. This integrated intensity, Lhk , is related to the modulus squared 
of the FF, |^^(q^fc)|2, by US] 

/,. = 2^y0(J)'-^|i^(q,,)r (1) 
\ 4 / <I>g qhk 

Here, V is the illuminated sample volume, (fi is the incident neutron flux density, 
A,i is the neutron wavelength, 7 is the magnetic moment of the neutron in nuclear 
magnetons, and $0 = h/2e is the flux quantum. General predictions from Ginzburg- 
Landau (GL) theory, whether by numerical methods or algebraic approximations, are 
for a monotonic decrease of the FLL FF with increasing field [27H30] . Although strictly 
valid only close to T^,, these models have been widely used and give a good description 
of results from many superconductors at lower temperatures j31H37| . Quasiclassical 
theory, which is valid well away from Tc, gives qualitatively similar results for both s- 
and d-wave orbitally-limited superconductors [35] • We will see that the conventional 
picture provided by all of these theories lies in stark contrast to the results that we 
now report. 

In figure 1 we show the field-dependence of the FLL FF for first-order reflections, 
at various temperatures between 50 and 1250 mK. At all temperatures, the FF initially 
rises with increasing field before reaching a maximum, and then begins to fall again on 
approaching Hc2- At temperatures of 500 mK and below, the FF remains finite all the 
way up to Hc2, where the FLL signal disappears abruptly upon entering the normal 
state. This is consistent with the first-order nature of the superconducting transition 
seen in thermodynamic measurements at low temperatures [II H], and predicted 
theoretically for strong Pauli- limiting [S]. Above 500 mK, the field-dependence of the 
FF becomes increasingly conventional, smoothly approaching zero at Hc2 , as expected 
for a second-order transition to the normal state. However, even at these higher 
temperatures, Pauli-paramagnetic effects remain important, leading to a maximum in 
the FF at intermediate fields. 

Figure 2 shows the regions within the superconducting phase diagram occupied 
by the various FLL structures [53]. For the hexagonal structures, two of the six Bragg 
spots in each of the two domains lie along (110) directions. In the high-field region. 
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wc find that the value of |F(qj^Q)|^ for these spots is ~ 20% larger than for any of the 
other four spots. This FF anisotropy is consistent with that expeeted to arise due to 
the slight distortion of the hexagonal structure away from the isotropic hexagon [53] . 
As a consequence of this, the magnitude of qio is always slightly less for the two (110) 
aligned Bragg spots than for the other four spots. To accommodate this variation in 
intensity, the FF values shown in figure 1 are an average over all six spots. Whenever 
the FLL structure changes, both the magnitudes and directions of q/jj, alter; however, 
figure 1 shows that any changes in the value of |-F(qio)P at the FLL structure phase 
boundaries are no bigger than the errors. Finally in figure 1, we see that at low field 
there is a common value of FF, independent of temperature, showing the differences 
between the data at different temperatures occur at higher fields. We may interpret 
this behavior as showing that the effects of Pauli limiting are small at low fields, 
and that for T < Tc/2, orbital-limiting effects are fairly temperature- independent. 
However, at high fields, Pauli-limiting is dominant and the proximity of Hc2(T) causes 
major changes in the behaviour of the FF. 

In figure 3 we show the results of investigating the detailed temperature- 
dependence of the FF in the high field hexagonal FLL structure phase. For the fields 
of 4.60 T, 4.85 T and 4.90 T, we have recorded the temperature-dependence of the 
diffracted intensity of a first-order Bragg reflection. At the two higher fields, these 
temperature scans pass through Hc2 almost parallel to the upper critical field phase 
boundary; nonetheless, there is a sharp fall in intensity as the first-order boundary 
is crossed, with perhaps a slight smearing arising from tiny differences between 
crystallites in the mosaic. We estimate that the contribution due to demagnetisation 
effects to the width of the discontinuity at Tc2 is unimportant. At the lower field, a 
slower variation with temperature of the normalised intensity is seen, reflecting the 
second-order transition at Hc2- However, as shown by the absolute data in the inset, 
the FF is at its maximum at the lower field, deeper within the mixed state, and falls 
as either field or temperature is increased towards Hc2- 

Previous discussion of the field-induced amplification of the FF emphasised the 
inability of conventional theory to explain our results at low temperatures |23] . Here, 
we consider the FF behavior in the context of an extension [55] to the quasi-classical 
Eilenberger theory [39]. In this work, Pauli paramagnetic effects are incorporated 
by adding a Zeeman energy term fiB{r) to the Eilenberger equations [25], where the 
parameter /i represents the relative strength of the paramagnetism |40j . This theory 
was first used to successfully describe the similar, but less extreme, field dependence 
of the FF in TmNi2B2C [JT]. At low temperatures, and for large values of ^, the 
model predicts not only an initial increase of the FF with increasing field, but also a 
decrease close to Hc2 [25]. To obtain a reasonable qualitative agreement between the 
theory and the observed FF field-dependence in CeCoIns requires /i ~ 2.6. However, 
experimentally determined material parameters for CeCoIns imply a ~ 1.7 for H || c. 
Hence, for a quantitative agreement with our data, some refinement of the theory is 
required. For instance, it might be necessary to incorporate a more three-dimensional 
character of the Fermi surface of CeCoIns [42], or account for the non- Fermi liquid 
behavior which is most prominent near Hc2- We further mention that heat capacity 
measurements [43] reveal a field-dependence to the low energy density of states (DOS), 
which further needs to be accounted for [33]. The inclusion of a field-dependent DOS 
would be consistent with the observation by thermal conductivity studies of two- 
band superconductivity in CeCoIns [45], which imply that some of the supercarriers 
become depaired at very low fields. These depaired carriers would not only be available 
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to contribute to the paramagnetic response at high fields, but their absence at low 
fields might also explain why the form factor at very low fields is not given by the 
expected expression (see [23p using the zero- field value of the London penetration 
depth, Al [TUIHS] . An independent bulk measurement of Al at finite fields would test 
this suggestion. 

It was previously suggested that the FF maximum at 50 mK could mark the 
onset of the FFLO state [23]. Fig. 2 shows the temperature dependence of the field 
at which a maximum in the FF is observed, superposed onto the FLL structure phase 
diagram. It is clear that the FF maximum does not follow the FFLO phase boundary 
proposed experimentally |14ll8j or predicted theoretically jT7ll48] . Furthermore the FF 
maximum does not appear to be correlated with the high-field rhombic to hexagonal 
FLL symmetry transition. To understand the observed FF behavior close to Hc2, we 
therefore return to the the theory of Ichioka and Machida [25l , which reproduces the 
observed drop of the FF without including any effects of an FFLO state. Beginning 
at low fields, they predict an increase in the magnitude and spatial extent of the 
paramagnetic moment in the flux line cores, arising from the spin-split quasiparticle 
states located there. The spatially varying paramagnetic moment adds to the orbital 
contribution to the field variation in the mixed state to give an increase in the field 
contrast and hence in -F(qio)- The size of the flux line cores, as measured by the extent 
of the region with a suppressed order parameter, also increases, showing the effects 
of paramagnetic depairing in these regions. However, at high fields where the width 
of the expanded flux line cores becomes comparable with the flux line spacing, the 
suppression of the order parameter allows the paramagnetic moment to spread further 
through the entire flux lattice unit cell. This leads to an additional paramagnetic 
moment of the whole crystal which is observed in magnetisation measurements [5] , but 
to a reduction in the field contrast, as measured by SANS, and hence to a reduction 
in F(q]^g). We emphasise that these effects are a consequence of paramagnetic effects 
dominating pairing in the near-core regions where the order parameter is suppressed, 
and that these effects spread out from there. Thus the decrease of the FF close to Hc2 
is quite different in origin from the behavior seen in conventional materials. 

Next we discuss the abrupt fall of the FF seen at Hc2 and low temperature. The 
finite value of the FF at Hc2 can be directly related to the jump in the magnetisation 
at the upper critical field [2]. First, we note that the spatial field modulation 
B{r) = -^(l/ifc) 6'*^'"'= where the sum is over all FLL reciprocal lattice indices 
h and k, and F(q = 0) is equal to the average induction {B). Close to Hc2 it is a 
good approximation to include in the sum only the dominant Fourier component 
-^(*l{io}); which is the FF measured in our SANS experiments. In the case of 
CeCoIns it is important to repeat that the field modulation contains contributions 
from both the orbital screening currents and the polarisation of the unpaired electron 
spins. However, just below Hc2j the spin magnetisation is by far the dominant 
contribution to the field modulation, such that B{r) = /io[Af(r) -I- -ff^ppl ] S5' 
Assuming that the magnetisation in the "normal" flux line cores is equal to -^normal 
just above H1.2, the average magnetisation in the mixed state is easily found to be 
(fioM) = ^^normal ~ 6|F(qjQ)|. This is illustrated in Fig. 4. Our measured value 
of the FF of 0.4 mT just below Hc2 at 50 mK would thus correspond to a jump in 
magnetisation of A{fioM) = 2.4 mT as the field is increased through the upper critical 
field. This is in very satisfactory agreement with the results of bulk magnetisation 
measurements, which give a jump at Hc2 of 3 mT [2]. 
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We now consider the effects of a flux line core paramagnetism on the sequence 
of FLL structures shown in Fig. 2, for fields H > 0.5Hc2 and low temperature. We 
note that this sequence of transitions from square to rhombic and then to high-field 
distorted hexagonal, is the exact reverse of what is observed for H < 0.5iJc2- In the 
central square phase, the FLL nearest neighbor directions are aligned with the nodes of 
the dx2-y2 gap; this is consistent with the low-field sequence of transitions being driven 
by flux line core anisotropy, which reflects that of the order parameter. However, on the 
approach to Hc2, paramagnetic depairing causes the flux line cores to begin to expand 
and also become more isotropic [H]. We expect that this suppresses the anisotropy 
which stabilised the square phase, and that core expansion is the driving mechanism for 
the high-field sequence of FLL structures. This may also explain the similar sequence 
of FLL transitions observed in TmNi2B2C [49], where Pauli paramagnetic effects 
have also been demonstrated [41]. At higher temperatures, although the high field 
effects of paramagnetic depairing on both the FLL FF and the FLL structure become 
suppressed, figure 2 shows that the re-entrant square to rhombic transition persists to 
temperatures above those for which Hc2 is first-order. This is a further demonstration 
that paramagnetic effects remain important at higher temperatures where Hc2{T) is 
becoming orbitahy-hmited. These deductions are consistent with those of another 
recent study based on a GL type theory, that also reproduces the observed sequence 
of high-field transitions and ascribes them primarily to the effects of paramagnetic de- 
pairing and Fermi surface anisotropy [5^. However, we point out that GL theory is 
not expected to be numericahy accurate over this region of the {H, T) phase diagram. 

4. Conclusion 

In conclusion, we have studied the detailed field- and temperature-dependence of the 
diffracted neutron intensity from the flux line lattice in CeCoIns. At high fields and 
low temperatures, we most clearly observe the effects of Pauli paramagnetism in the 
flux line cores, which results in an increase of the FLL form factor that extends well 
into the temperature region where Hc2 becomes second order. Close to Hc2 the flux 
line lattice form factor decreases, which we attribute to the effects of a paramagnetic 
suppression of the Cooper pairing. We suggest that further consequences of this 
paramagnetic de-pairing are manifested in the observed high field sequence of FLL 
structure transitions. 
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Figure 1. Field-dependence of the square of tiie fiux line lattice form factor 
for first-order Bragg reflections, li^(qio)Pi for temperatures up to 1250 mK. The 
dashed lines are guides-to-the-eye. The points at which the form factor goes to 
zero are obtained from measurements of Hc2 I14| . 
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Figure 2. Superconducting piiase diagram for CeCoIns with H ||c, indicating 
the regions occupied by the different fiux line lattice (FLL) structures reported in 
I2m24| . Phase boundaries separating regions of different FLL structure are shown 
by gray dotted Unes. The inset diagrams show the three FLL configurations. Here 
the crystaUinc [100] direction is vertical, and first-order Bragg spots belonging 
to different FLL domains are denoted by open and filled circles respectively. 
The green and black solid lines show where the superconducting to normal state 
transition is first- and second-order respectively. The open blue circles show the 
temperature evolution of the maximum of the form factor (|i^(qio)P) , with the 
blue dashed line being a guide-to-the-eye. The error bars of the data points are 
comparable in size to the data symbol. The area enclosed by the green solid and 
red dash-dot lines provide an estimate of the region occupied by the FFLO phase, 
as deduced from data shown in |14II18| . 
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Figure 3. The normalised temperature-dependence of the diffracted intensity at 
the peak of the rocking curve for a first-order hexagonal flux line lattice Bragg 
spot lying along (110). Data were collected at fields of 4.60 T, 4.85 T and 4.90 T. 
The dashed lines are guides to the eye. The inset shows the same data as the 
main figure, but in absolute (un-normalised) units of counts per standard monitor. 
Since the rocking curve width is resolution limited, these are oc Iio in Eqn. (1). 
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Figure 4. A schematic of tiie spatial variation of tiie local magnetisation 
along the nearest neighbor direction between adjacent flux line cores, and for 
an applied field (i^appi ) just below Hc2- The magnetisation is expected to be 
close to the the normal state value in the cores, and decreases between the flux 
lines where the Pauli paramagnetism is suppressed by the Cooper pairing of the 
supercarriers. The local magnetisation varies from +6l-F(qjQ)| to — 2|i<'(qj^g)| 
relative to iio{M). The minimum in the local magnetisation is ^q{M) — 3|-F(qj^Q) 
and is obtained along the FLL next nearest ncightbor direction (unit cell diagonal) 
(not shown). The shading indicates the origin of the reduced magnetisation in 
the superconducting state. 



